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Abstract 

String-inspired cosmologies, whereby a non-singular curvature bounce is induced by a 
general-covariant, T-duality- invariant, non-local dilaton potential, are used to study 
numerically how inhomogeneities evolve and to compare the outcome with analytic 
expressions obtained through different matching conditions across the bounce. Good 
agreement is found if continuity across the bounce is assumed to hold for TZ, the curva- 
ture perturbation on comoving hypersurfaces, rather than for the Bardeen potential. 
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1 Introduction 



The idea that string theory, with its fundamental length scale, could resolve the big bang 
singularity by effectively limiting space-time curvature has led, during the last decade, to 
a variety of models [1] where the big bang represents a turning (rather than an end)-point 
in the history of the early Universe. All these models share the property of exhibiting a 
bounce in curvature, since the pre (resp. post)-bang eras are characterized by growing (resp. 
decreasing) space-time curvature. Hereafter we shall refer to this generic class of models as 
pre-big bang (PBB) cosmologies. 

Technically speaking, the PBB phase can replace the inflationary phase of standard 
slow-roll inflation since it allows the cosmologically interesting scales we observe today to 
start well inside the Hubble radius, and thus to have time to homogenize through causal 
microphysics. An advantage of PBB cosmologies [2] is that they are less sensitive to the 
so-called trans-Planckian problem since the proper frequency of all relevant quantum fluctu- 
ation starts very much below the Planck (but not necessarily the string) scale; the opposite 
is true in the case of slow-roll inflation, where all relevant proper frequencies are initially 
much above the Planck scale, formally causing linear perturbation theory to break down. 
Another difference is that, while in ordinary inflationary models string (or Planck)-scale 
physics is washed out by inflation (except for the above-mentioned possible imprint on ini- 
tial conditions) , the opposite is the case in PBB cosmologies since the high-curvature phase, 
the bounce, comes there at the end of inflation. 

If this latter property of PBB cosmology is appealing -in that it opens an observational 
window on short-distance physics- it also makes it mandatory to understand how the bounce 
itself may (or may not) affect different observables. Unfortunately, this is not an easy 
problem to solve since the details of the poorly understood mechanism responsible for the 
curvature bounce may afl'ect, in a crucial way, how different perturbations behave across 
the bounce itself. As an example, while the original investigations [3] argued in favour of a 
strongly tilted (blue) spectrum of adiabatic density perturbations, more recently arguments 
have been given [4] for much flatter (and possibly even scale-invariant) spectra. In toy 
models where the bounce is schematized as a discontinuity across a space-like hypersurface 
those diflFerent predictions can be shown [5] to correspond to diflFerent ways of implementing 
Israel's junction conditions. This has led to much debate in the literature lately (see for 
instance [6]). 

An alternative approach is to use, as a theoretical laboratory, a model of background 
which is everywhere smooth and regular, so that we can follow both analytically and nu- 
merically the evolution of the different fluctuations throughout. One can then replace the 
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smooth evolution by a sudden transition, and check how different prescriptions match with 
the "exact" results. We will work, in particular, in the context of a perturbative approach 
based on the string effective action, where it is known that curvature singularities can be 
eliminated either at large curvature/coupling by a combination of higher-derivative and loop 
corrections [7] or, in the low-curvature regime, by the contribution of a non-local dilaton 
potential [8, 9, 10]. In this paper we will concentrate on this second possibility, and we will 
consider regular cosmological models obtained by using a non-local potential put in by hand. 
Such a potential may arise, however, from dilatonic loop corrections in manifolds with com- 
pact spatial sections. The corresponding non-local potential energy provides a dominant 
correction to the regularization of the background when an effective dimcnsionally-rcduccd 
coupling grows large, and leads to bouncing solutions smoothly evolving (in the Einstein 
frame) from accelerated contraction to decelerated expansion. 

This paper extends the results of a previous letter [11] by investigating the frame- 
independence of the results, by performing the above-mentioned comparison with different 
matching conditions, and by including a perfect fluid source both in the background and 
the perturbation equations. However, we shall limit our detailed discussion of scalar pertur- 
bations to a particular example of gravi-dilaton background without matter sources. The 
discussion of more general examples is postponed to further work [12]. 

Our analytical and numerical computations suggests that the scalar spectrum associated 
to a non-singular background can be correctly obtained through an (approximate) matching 
procedure, provided the matching prescriptions are such as to guarantee a smooth crossing 
over the bounce of the variable TZ, representing spatial curvature perturbations on comoving 
(or, at large scale, on constant energy density) hypersurfaces, and not of the variable 
the gauge invariant Bardeeen potential. It is possible, however, that such a resTilt is due to 
the fact that the bounce is triggered by a potential V that depends mostly on the value of 
the dilaton: as such, it may not apply when the matching is performed through a singular 
bounce, as it seems to be the case in the ekpyrotic scenario [13]. 

The rest of the paper is organized as follows. In Section II we introduce the low- 
energy effective action to be used throughout this paper and discuss, in particular, the 
possible origin of the non-local dilaton potential. In Section III we derive, both in the 
string and in the Einstein frame, the (non-local) field equations that follow from the above 
action, and present classes of exact solutions describing regular homogeneous bouncing 
backgrounds. The last part of Section III, containing the derivation of the exact solutions 
with fluid sources, is not essential for the understanding of the subsequent sections, and 
may be skipped by those readers who are only interested in the discussion of the matching 
prescriptions. In Section IV we introduce the full set of scalar perturbation equations in the 
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Einstein frame, including the perturbations of the fluid sources, both in the uniform-dilaton 
gauge and in full gauge-invariant form. We also give explicit relations to the String-frame 
perturbation variables, for easier comparison with previous papers and results. In Section 
V we estimate analytically, and compute numerically, the scalar perturbation spectrum 
(inside and outside the horizon) for a typical example of regularized background, and we 
discuss the possibility of reproducing the same spectrum through an appropriate matching 
procedure across a sudden bounce. Section VI is devoted to some general comments and 
concluding remarks. 

2 Non-local dilaton potential and low-energy effective action 

It is well known that implementing a curvature bounce in cosmology must involve some 
unconventional equation of state for the dominant source of energy. Such a possibility 
can arise in principle from higher-order corrections in a', string-loop corrections, or both. 
Examples of such "exit" mechanisms have been widely discussed in the literature [7]. In 
this paper wc wish to implement a low-curvature bounce and therefore we will not appeal 
to the former kind of corrections. We shall argue, instead, that string-loop corrections, in 
the presence of space compactification, can provide an efficient and appealing regularization 
mechanism. 

In a D-dimensional manifold, if n < D — 1 spatial dimensions are compactified, at ener- 
gies below the compactification scale the effective lower-dimensional dynamics is controlled 
indeed by the "reduced" coupling: 

gled = al/Vn , (2.1) 

where = is the string coupling of the full D-dimensional theory (0 being the dilaton) 
and Vn is the n-dimensonal compactification volume in string units. In standard superstring 
compactifications from D = 10toZ) = 4, the reduced coupling gl = q^/Vq is nothing but 
the four-dimensional gauge coupling that controls the dynamics of the "zero-modes" with 
respect to the six-dimensional internal space, i.e. of the light four-dimensional degrees of 
freedom. As an example, it is known that the running of g^^ is determined by 5(4 itself 
through the infrared behaviour of zero-mode loops. It is also known that T-duality (i.e. 
large-small volume duality) is preserved by loops, provided the reduced coupling is kept 
fixed. 

In a cosmological set-up where the whole of space is compact, the zero modes become 
the homogeneous components of the various fields, and their self-interactions should be 
controlled by the D = 1 effective coupling gl = gl /Vg. In analogy with the previous example 



3 



we may expect loop-corrections that preserve T-duality where, this time, e'^^*'' /Vg = e'^^^^ = 
gf is kept fixed. In particular, if a non-trivial effective potential is induced, it should satisfy 
such a symmetry and must therefore be of the form K// = V{Tp). 

It is known [8, 9, 10] that such potentials can be very useful for inducing curvature 
bounces. However, for a long time, it was not clear how to formulate such a framework 
in a fully general-covariant way. This problem was solved in [11], where the following 
generalization of e"'^^^^ to a function of all the coordinates x was proposed {d = D — l): 



Such a quantity, even if non-local, transforms as a scalar function of x under general coor- 
dinate transformations, and is automatically 0{d,d) invariant for a string-frame (S-frame) 
metric with d abelian isometrics. It should be remarked, however, that such a definition of 
Tp is only consistent for backgrounds such that d^ipd'^ip > 0. This will be always the case 
for all background solutions we will consider in this paper. 

In the spirit of the previous discussion, our toy model will consist of taking the tree-level, 
low-energy effective action of (say heterotic) string theory, and of adding to it a suitable 
effective potential depending only upon Tp(x) , while other possible matter fields are included 
into the effective action Sm- Our effective action will then have the form, in d+1 dimensions, 
and in the S-frame: 



^ = -7;^ I d''^^^\f\9\ \R + (V^')' + V{e-^)\ + S^. (2.3) 

Here Ag = = (27ra')^/^ is the string length scale, and the non-local function ip = '(p{(p) 

is defined in Eq. (2.2). Our conventions are: diag g^i^ = (-1 . . .), R^ua ^ = d^Ti^a ^ + 

r I3r p — and E = E ^ 



3 General field equations and background solutions 
3.1 Field equations in the string Frame 

Starting from the action (2.3) the field equations are obtained, as usual, by taking functional 
derivatives of the action with respect to the metric and the dilaton. By adopting the 
standard definition of the matter stress-energy tensor, 

oVifl ^Mi" (3-1) 



Sgt'^ix) 2 
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we note that 

5 



6g^^•'{x) 



where 



and 



V' = 



dV _ ^dV 



(3.3) 



(3.4) 



We are using the convenient notation in which an index appended to round brackets, (. . .)j;, 
means that all quantities inside the brackets are functions of the appended variable. Simi- 
larly, = ip{x). 

Collecting all terms, and multiplying by e"^ / vl^Ii we obtain the final equation 

G^r. + V^V^^ + [iy^f - 2V V - V 
where Gnv is the usual Einstein tensor, and 



Je-^^(a^ 7^.^i = Af-^e^T^., (3.5) 



h{x) = d'^+'y F')^ 5 (^, - ^y) . 



(3.6) 



We now take the functional derivative of the action with respect to the dilaton, and 
define the scalar charge-density of the matter sources, a: 

6Sin 



(3.7) 



The differentiation of the non-local potential gives 
S 



j J d^+^x' (^e-^y) ^ 6{^, - ^,,), (3.8) 



5 



where 5' denotes the derivation of the delta distribution with respect to its argument. The 
last term of the above equation can be rewritten as 



-2ip 



while the term containing the 5' integral can be rewritten as 

'dV 



where 



h{x) = d'+\ 7(0^) 5' - ^y) 



(3.9) 



(3.10) 



(3.11) 



Collecting all terms, and multiplying by e'^/ \/\g\, we finally obtain the dilaton equation to 
read: 

R + 2VV - {V^f + V-^ + e-^^^/i - e-'^V'h = Xt^e^a, (3.12) 

dip \/{oipY 



where we have defined 



(3.13) 



Eqs.(3.5) and (3.12) were already given (without presenting their explicit derivation) in 
[11]. In the absence of matter sources they are not independent, being related by 



5S 



1^ f^S 



(3.14) 



For later applications, it may be useful to consider also a combination of the metric and 
dilaton equations, by eliminating the scalar-curvature term appearing in Eq. (3.5). The 
result is 



h = 



(3.15) 



Let us finally stress the difference with the more conventional case in which the potential 
appearing in the action is a local function of the dilaton, V = V{(p). In that case, the vari- 
ation of the action provides a set of scalar-tensor equations which can be directly obtained 
from Eqs. (3.5), (3.12), simply by neglecting the terms with the Ii and I2 integrals, and 
replacing in Eq. (3.12) dV/dtp by dV/dip. 
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3.2 Field equations in the Einstein Frame 



In this paper we will complement the analysis of [11] by studying the evolution of metric 
perturbations directly in the Einstein frame (E- frame), where the dilaton is minimally 
coupled to the metric and canonically normalized. To this purpose, it is convenient to write 
down the field equations starting directly from the E-frame action. This is obtained from 
the S-framc action (2.3) through the conformal rescaling (here we will denote with a tilde 
the E-frame variables): 



9nu 



2^/{d-l) 



(3.16) 



We have explicitly introduced the Planck length scale Ap to restore the conventional canon- 
ical normalization of the gravitational variables. Prom the action (2.3) we thus obtain, up 
to a total derivative, 



2A' 



d-l 



where Ap defines the D = {d + l)-dimensional gravitational coupling Gd as 

1 



and where 



By introducing the (E-frame) canonically normalized scalar field ip: 



(3.17) 

(3.18) 
(3.19) 



Ml 



ip = ^ 



the action can be written in canonical form as 



-1 



R 



1 



IGttGd 2 



\ 2 



V 



+ Smig,'p), 



where 



V 



(3.20) 



(3.21) 



(3.22) 



We shall work in units I&'kGd = 1, i-e. Mp ^ = 2. In this case eq. (3.20) reduces to 

= ap, c=\ . (3.23) 



d-l 
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We shall also absorb the factor (A^/Ap)^ into the function F(^), with the understanding 
that now the dilaton potential is rescaled from string units to Planck units. The action 
finally becomes 



S{g,^) = j d<'+^x^\ [-R+\ (V^)' - V 



(3.24) 



where 



V{g,(p) = e'^V{e-n, 



In order to obtain the field equations we will now take functional derivatives of the above 
action with respect to g^^^ and Ip. For the matter action Sm we set 



5S„ 



(3.26) 



and find that, according to the definition of g and the E-frame sources terms are related 
to the S-frame ones by the rescaling 



2 c 



(3.27) 



By separating explicitly the non-local part V of the potential V of Eq. (3.25), and following 
the same procedure as in the S-frame, the variation of the full action with respect to g^'^ 
gives 



1 ~ 



1 



1 



1 



(3.28) 



where 



Ji{x) = ^J d^+^y y\^\V'^^5{^. - (?,), (3.29) 



and, as before, V = {dV/de~'^). The variation with respect to ip gives 



V'^ip +-a + -cT + e'"^ 



_dv 

2^ 



W^Jl - V'J2 



where 



= 0, (3.30) 



(3.31) 
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The presence of the trace of the stress tensor, in this equation, is due to the y?-dependence 
of the matter action through the rescaled metric g, leading to the variational contribution 

The combination of the two equations (3.28), (3.30) leads to the generalized conservation 
equations 

V.f;, ''-^(a + cf) V^^ = 0. (3.33) 

We finally recall, for comparison, that the (more conventional) equations obtained from 
the Einstein action (3.24) in which the S-frame potential V is a local function of the dilaton, 
V = V{ip), can be simply obtained from Eqs. (3.28), (3.30) by dropping the integrated 
contributions Ji, J2, and by replacing —{c/2){dV/d(p) by dV/dip in Eq. (3.30). 



3.3 Homogeneous background solutions 

For the cosmological applications of this paper we have to start with (and then perturb) a 
homogeneous, isotropic and spatially flat background. Considering first the S-frame equa- 
tions we set: 

g^jy = diag(l, -a^6ij), " = diag(p, -p6ij), tp = ip{t) (3.34) 

(we are assuming that the matter stress tensor can be expressed in the perfect fluid form), 
and we obtain we obtain, from Eq. (2.2), 

g-^ = e-'^a'^, (3.35) 

where, assuming spatial sections of finite volume, we have normalized the scale factor a as 
to absorb in a"^ the S-frame spatial volume in string units. 

Let us write the associated S-frame field equations, in units 2Af~^ = 1. The (0,0) 
component of eq. (3.5) gives 

1^ -dH^-V = e^p, (3.36) 
while the (ij) component of eq. (3.15) gives 

H-Ht=l-e'^(p-la], (3.37) 



2 V 2 

and the dilaton equation (3.12) gives 



2ip - ^"^ - dH^ + V - ^ = le'^a. (3.38) 

olp 2 
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Their combination leads to the matter conservation equation 



'p + dHp=ha{ip + dH) , (3.39) 



where p = pa'^, p = pa'^, a = cra'^. 

When written exphcitly in terms of (p, the above equations can be easily compared with 
the homogeneous cosmological equations associated with a local potential V{ip). Besides 
the replacement dV/dlp — > dV/dtp in the equation for the dilaton field, there is a second 
remarkable difference appearing in the spatial equation corresponding to Eq. (3.37). For 
a local potential V{ip) such equation contains, on its left-hand side, the additional term 
{l/2){dV/dip). As illustrated below, this difference is responsible for making our equations 
admit exactly integrable and regular bouncing solutions. 

In the absence of matter sources, T^j^ = 0, a = 0, the dilaton equation (3.38) becomes 
a consequence of Eqs. (3.36), (3.37) (provided Tp ^ 0). The whole system of background 
equations can then be reduced to quadratures, namely 

H = me^ t = p dX' [d + A'V(mA')] , (3.40) 

where m is an integration constant, and we have defined A = m''^e~'^. As noticed in [8], it 
is then possible to obtain, quite generically, non-singular solutions if the function appearing 
between square brackets in the above integral has a simple zero. Furthermore, if VX'^ —>■ 
at large |A| (corresponding to a potential generated beyond the two-loop level), then the 
solutions asymptotically approach, at large \t\, the usual isotropic solutions dominated by 
the dilaton kinetic energy (see for instance [1]), characterized by H\t\ = 1/y/d. 

To illustrate this point, and to provide explicit analytical examples of regular back- 
grounds, we shall consider here the particular class of potential defined by 

V{mX) = A-2 f fa - (mA)-H ^'^^"^ - d) , (3.41) 



and parametrized by the dimensionless coefficients a and by the "loop-counting" parameter 
n. For Q > and ra > 0, Eq. (3.40) can then be integrated exactly, and leads to the general 
solution 

H = me'^ = m 



a 



l/2n 

.l + (amt)2". ' ^^-^^^ 
whose "bell-like" shape describes a bouncing evolution of the curvature scale, i.e. a smooth 
transition from growing to decreasing curvature, and from accelerated to decelerated ex- 
pansion. 

The asymptotic behaviour of this class of solutions is controlled by the parameter a. As 
\t\ —>■ GO one finds indeed H = a(^~^")/^"/|t|, so that the background (3.42) evolves from 
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an initial accelerated, growing dilaton configuration, parametrized for t — > — oo as 



a ~ {-ty 



{l-2n)/2n 



ip rsj — 



d a(l-2")/2n + 1 



ln{-t), 



t < 0, 



to a final "dual" and time-reversed configuration, parametrized for t — ^ +00 as 



,a{l-2n)/2n 



t > 0. 



(3.43) 



(3.44) 



The minimal, pre- and post-big bang solutions dominated by the dilaton kinetic energy [1] 
are thus recovered, asymptotically, for 



a 



^n/(2n-l)^ 



(3.45) 



If this condition is not satisfied the background is potential-dominated even asynptotically. 

For such a class of bouncing solutions the dilaton keeps growing ast —>■ +00, if a^^'^'^^^^'^d 
> 1. It is thus possible that a solution, which is regular in the S-frame metric corresponding 
to the action (2.3), acquires a curvature singularity when transformed to the E-frame metric 
corresponding to the action (3.21). Let us compute, indeed, the E-frame curvature scale, 
or the Hubble parameter 

din a 



H 



dt 



H 



d-1 



dt = dt e-'^'^'^-^\ 



(3.46) 



where t is the cosmic-time parameter in the syncronous E-frame gauge, while the dot still 
denotes differentiation with respect to the S-frame cosmic time parameter. Asymptotically, 
at large positive times, we obtain from Eq. (3.44): 

Jl-2n)/2n" 



H 



d- 



^_4_[„(l-2n)/2„_i]^ 



(3.47) 



The solution is thus bouncing even in the Einstein frame (irrespective of the possible growth 
of the dilaton), provided 

< 1. (3.48) 



a 



(l-2n)/2n 



This conditions guarantees that the Hubble factor of the E-frame metric approaches zero 
from positive values as f — ^ -|-oo, and we note, in particular, that the condition is always 
satisfied by the free vacuum solutions defined by Eq. (3.45). 

It is important to stress, for subsequent applications to metric perturbations, that the 
asymptotic form of the Einstein scale factor a, when written in conformal time, becomes 
completely independent of the parameters a and n: 



a ^ ae 



a ^ ae 



i-v) 



1 

d-1 



-f/{d-l) 



r? < 0, 
r?>0. 



(3.49) 
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where drj = dt/a = dt/a. The two asymptotic branches of 0(77) are thus simply related 
by a time reversal transformation, as appropriate to a "self-dual" background [14] in the 
Einstein frame. The Einstein scale-factor 0(77), on the other hand, represents the "pump 
field" determining asymptotically the spectrum of scalar and tensor metric perturbations 
[15] (as also illustrated in the subsequent sections). All the solutions of the considered class 
are thus expected to lead to the same asymptotic spectrum of amplified metric fluctuations, 
irrespectively of the values of a and n. For the analytical and numerical computations of 
this paper it will thus be enough to use the particularly simple case n = 1, a = d, m? = Vq 
of the general class (3.42), corresponding to the (non-local) two-loop potential V = —Vqc^'^. 
In such a case the exact solution, for d = 3 spatial dimensions, can be explicitly written as 
[10] 

1/V3 I 



air) 



+ Vr2 + 1 



T = t/to, t^^ = e^^^^/Vo, (3.50) 

where (without loss of generality) we have set a = 1 at the bounce, occurring at t = 0, 
and where (/jq is an integration constant. Such solution is characterized by the two relevant 
parameters, Iq^ and e'^^ , corresponding (up to numerical factors) to the Hubble parameter 
and string coupling at the bounce. 

For later applications it is now convenient to present also the homogeneous limit of the 
E-frame equations (3.28), (3.30). Considering a homogeneous, isotropic and spatially flat 
E-frame background we obtain, from Eq. (3.19) and (3.29), 

dV 



where we have absorbed, as before, the total, finite volume (now in Planck units) of the 
compact spatial sections. The (00), {ij) components of the Einstein equations (3.28), and 
the dilaton equation (3.30), thus become, respectively (in units 2Ap~^ = 1): 

d(d-l)#2 = p + l(^%e^^y, 

2{d - \)H + d{d - 1)#2 = -p- - e^?y) - 6^=?^ 



2^ ; dip' 

-•1 c -/ cdV\ 

(p + dH(p+-a+-{p-dp)+ e''^ [cV - -— j = Q, (3.52) 

where the dot now denotes differentiation with respect to the E- frame cosmic time t. One 
can easily check that such equations can also be directly obtained from the corresponding 
S-frame equations (3.36)-(3.38), by performing the canonical rescaling 

a^a = ae-'^/('^-^) , dt ^ dt = dte"^'^'^'^^ , ip^^ = c^, 

~ d+l ^ d+1 ^ (J d+1 

p ^ p = pe'^-^^ , p —>■ p = ped--^"^ , a —>■ a = —ed--^'^ , (3.53) 
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connecting the two frames. The combination of the three equations (3.52), in particular, 
leads to the conservation equation 

P+dH{p + p)-^{p-dp)^-^a^ = 0, (3.54) 
and to the useful identity: 

~ -2 dV ~ 

2{d-l)H + p + p + if +^e"'^ = 0. (3.55) 

(j(p 

It is important to compare these equations with those obtained from an action identical 
to (3.24), but with a local potential V = e'^'^V{(p). The (00) equation and the conservation 
equation are the same, but the other equations are different. One obtains, for a local 
(S-frame) potential V{(p), 

d(d-l)#2 = p+l~%e^?y^ 
2{d - l)H + d{d - 1)H^ = -p- - e^^V^ , 

^ + dH(p + -{p-dp) + e^'^ [cV + — ] = 0. (3.56) 



2 2'"^ V dip 

The combination of the first two equations, in particular, leads to the usual identity 

2{d-l)H + p + p + ^'^ = 0, (3.57) 

SO that the behaviour of H is completely determined by ip and by the fluid variables. For 
the background solutions considered in this paper this is no longer true, because of the 
explicit contribution of the potential to Eq. (3.55). 

The (E-frame) background equations (3.52) will be applied, in the following sections, 
to discuss the evolution of scalar metric perturbations, and to compute their final spectral 
distribution. In d = 3 spatial dimensions, and in terms of the frame- invariant conformal 
time parameter rj, they can be written explicitly as 

67^2 = pd^ + ^(^'2 + e^F5^ (3.58) 

+ 2^2 = -pd^ - Q<^'2 _ g?^~2^ _ e^^d^ (3.59) 

^" + 2W + ^aa'' + l{p-Sp)d^ + e^(v- d^ = 0, (3.60) 

where 7i = d' /d, and the prime denotes differentiation with respect to rj. The conservation 
equation and the identity (3.55) read, respectively, 

p' + 3n{p + p) - \ip'{p-m - = 0, (3.61) 
iin' - n^) + + (p + ^)a2 + e?^S2 ^ Q_ (3_g2) 

CHp 
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It is finally worth noticing, even if not essential for the understanding of our subsequent 
discussion on the matching of perturbations, that regular bouncing solutions are possible 
even in the presence of matter sources. Working with the S-frame equations we may con- 
sider, in particular, the special case in which dV/d(p = 2V, and in which the homogeneous 
equations (3.36)-(3.38) can be reduced to quadratures for a wide class of equations of state 
satisfying the condition that the ratio p/p is an integrable function of the time-coordinate 
X, related to cosmic time by [9] 

2dx = Lpdt, L = const. (3.63) 

We can take, for instance, 

V = -V^e^, P = 7P, 0" = JqP, (3.64) 

with Vb, 7, 7o constants. In such a case, assuming that 7 7^ 70/2 and that 70/2 7^ 1, a first 
integration and a combination of the above equations leads to [16] 

-/ 0/^1 lo\x + xo a' ( -iQ\x + xi 

where 

D{x) = L^e-^p = (^1 - (x + xof -d(^j-^^\x + xif + LVq, (3.66) 

and xi,xo are integration constants. A second integration may then provide regular solu- 
tions for a and Tp if D(x) has no real zeros, namely if the discriminant of the quadratic 
form D{x) satisfies < 0. This is possible provided Vq > 0, i.e. for V{'ip) < 0. In the 
simple case 70 = we get, for instance, 

A^ = d-f\xi - xof + L^Voidj^ - 1), (3.67) 

so that the condition A^ < can be implemented for a set of reasonable equations of state, 
characterized by dj^ < 1. 

In the special case a = and p = we are lead, with a first integration, not to eqs. 
(3.65), but to the following set of equations [16], 

^ = _2^±^, ^ = ^, D = {x + xof-dxl + L'Vo. (3.68) 

Regular solutions are thus allowed for A^ = dx^ — L?Vq < 0. The choice xi = 0, Xq = 0, 
leads to the particularly simple exact solution describing a trivially flat (String-frame) space- 
time, with the dilaton performing a (time-symmetric) bouncing evolution sustained by the 
presence of a constant energy density of the (dust) matter sources: 

a = ao, P = Po, P = 0, = /■/, no ' (3-69) 

1 + [t/to) 
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where oq, po are integration constants, and 



e 



4 



/3o = T2 =^oe'^°- (3-70) 



/2 



It should be noted that such solution acquires a less trivial representation in the Einstein 
frame where, by using the appropriate transformations (3.53) one finds, in conformal time. 





where rjo = to/ao. In this frame the metric smoothly evolves from accelerated contraction 
to decelerated expansion, sustained by the dilaton and by a non-trivial, "bell-like" evolution 
of the energy density of the pressureless fluid. The two asymptotic branches of such a "self- 
dual" solution are simply related, in the Einstein frame, by a time-reversal transformation, 
like in the previous class of vacuum solutions (3.49). This simple background, together with 
a more general class of solutions, will be used in a forthcoming paper [12] to discuss the 
spectrum of metric perturbations in the presence of perfect fluid sources. 



4 Perturbation equations in the general case 
4.1 Perturbation equations in the Einstein frame 

In this section we will perturb the E- frame background equations presented in subsection 3.2, 
by expanding the background fields around the homogeneous configurations of subsection 
3.3, in order to obtain the linearized equations governing the evolution of the (first order) 
gravi-dilaton and matter fluctuations. Note that, from now on, we will drop the tilde for 
the sake of a simpler notation, but all variables appearing in the subsequent equations are 
to be interpreted (unless otherwise stated) as E-frame variables. 

The computation of the tensor (transverse and traceless) part of the metric perturbations 
can be easily performed in the syncronous gaiige, by setting 

500 = 1, 50^ = 0, gij = -a^5ij, 

Sgoo = = 6gio, dgij = hij, g'^hij = = djhi^, (4.1) 

while keeping the dilaton and the matter sources unperturbed, Stp = 0, ST^i^ = = da. We 
then find that all the perturbed terms of Eqs. (3.28), (3.30) are vanishing, except 

SRi ^ = -^Ohi^ = (k'+dHhi^ - . (4.2) 
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For each polarization mode h we thus recover the usual (E-frame) evolution equation, which 
can be written in Fourier space (using conformal time and in cZ = 3) as 

h'l + 2neh'k + k'^hk = 0. (4.3) 

The associated canonical variable [15] 

^^ = ^ah, - (V^ + ^) = 0, (4.4) 

diagonalizcs the kinetic term in the effective action for tensor metric perturbations, has 
canonical dimensions and can be correctly normalized to an initial spectrum of quantum 
vacuum fluctuations {\nk\ |2fe|~^/^ for r] — oo). For a thorough discussion of tensor 
perturbations see [11]. 

To determine the evolution of scalar perturbations in a conformally flat background, 
9ij,u = diag a^(l,— 1,— 1,— 1,...), we will keep using the conformal time coordinate, and we 
will work in the so-called "uniform-dilaton" gauge [11], by setting 

S(p = 0, SgoQ = 2a^(t), Sgij = 2a^(V'(5y - didjE), 

Sn'' = Sp, STi^ = -Sp6{, STi'' = {p + p)diu, (4.5) 

where u is the usual velocity potential for the fluid perturbations [15]. We find, in this 
gauge, 

,57;: = 0, (5(e-^) = 0, <5/i=0. (4.6) 

The perturbation of the Einstein equations (3.28) leads then to 

SG", = 1st;: + \5g^^d,^d^^ - (Sg'^^d^^dn^) + \^''ll^4>. (4.7) 

where the last term of this equation has been computed using already the homogeneity 
assumption = (/^(ry). 

Let us consider the above equation for the case d = 3. The [i ^ j) component gives 

E" + 2nE' + V - <^ = 0. (4.8) 

The (zO) component gives 

i^' + n(l>=^ip + p)a\. (4.9) 

The (00) component gives 

+ HE') - m{ip' + n(i>) = ^a^Sp - (410) 
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The (i = j) component gives 

{2n' + n'^)^ + H<t>' + V" + W-iP' = ^a^5p - ^(^'2</) - ^e^^aV (4.11) 

(we have also used the (i 7^ j) condition, i.e Eq. (4.8)). 

Consider now the Hmit in which the fluid sources (and fluid perturbations) are absent 
(or negligible). In that case, the first two equations are reduced to 

E" + 2nE' + ^-(l) = 0, (4.12) 
i;' + H^ = 0, (4.13) 

Eq. (4.10), using the momentum constraint (4.13) and the background equation (3.58), can 
be rewritten as 

2V2(V' + HE') - eUip' - e^Va^cf) = 0. (4.14) 

Eq. (4.11), using the background equation (3.58) for (^'^, and the background equation 
(3.59) for dV/dlp, can be rewritten as 

V'" + 2Hil)' + n4>' + {W + 2K^)(/) = 0. (4.15) 

It can be easily checked that the four above equations exactly correspond to the E-frame 
transformed version of the perturbation equations computed in the S-frame in the absence 
of fluid sources, already presented in our previous paper [11]. 

In the presence of fluid sources, the system of equation (4.8)-(4.11) can be completed 
by perturbing the matter conservation equation (3.33). For the application of this paper 
we will limit ourselves to the case a = 0, c^ct = 0. The time component of the perturbed 
equation gives then 

5p' - (p + p)V'^u +{p + p){V'^E' - 3^-') + m{5p + 5p) = ^<4>'{5p - 36p), (4.16) 
while the space component gives 

u' + mu+^-^u-(f>=-^. (4.17) 
p + p p + p 

4.2 Gauge- invariant scalar perturbations 

The perturbation equations written in the previous subsection have been obtained using 
the particular gauge (4.5), and are not invariant under the infinitesimal coordinate trans- 
formations 

ri^ri + e^{r],x'), x' ^ + d'e{r],x'), (4.18) 
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still preserving the scalar nature of the metric fluctuations. However, even within the more 
general decomposition of dilaton and scalar metric perturbations, i.e. 



Sgoo = 2aV, Sgij = 2a^{tljSij - didjE), (4.19) 

and for any given choice of gauge for the variables {x,(pi '4'i B,5p,5p,u}, it is always 
possible to define an independent set of gauge- invariant perturbations. For instance [15], 
the set of variables X, £, P, W} defined by 

^ = (I) + - [{B - E')a\' , ^ = '^-H{B-E'), 

X = x + ^'{B-E'), £ = 6p + p'{B-E'), 

P = Sp + p'{B-E'), W = u + B-E', (4.20) 

is invariant under the infinitesimal transformation (4.18). The gauge invariant Bardeen 
potentials ^ and in particular, coincide with the scalar components of the metric fluc- 
tuations in the so-called longitudinal gauge E = = B [15]. We also introduce, for later 
use, another useful gauge invariant combination, 

n = -4^-^x = -^- ^x, (4.21) 

representing (in the absence of fluid sources) the perturbation of the spatial part of the 
scalar curvature [17] with respect to comoving hypersurfaces, i.e., spacelike 3-d surfaces 
characterized by a uniform dilaton distribution. In our case, the perturbation variables 
associated to the uniform-dilaton gauge can be expressed in terms of the variables (4.20) as 
follows: 

, X n n X' ip" 

(fr (p' (p' (p' (pr 

dp = £-^X, Sp = P-^X, u = W-^. (4.22) 
(p' ip' (p' 

By inserting these deflntions into the previous set of equations (4.8)-(4.11), (4.16)-(4.17), we 
can thus obtain a linearized set of perturbation equations written in a fully gauge-invariant 
form. In particular, Eq. (4.8) becomes 

$ = *. (4.23) 

Eq. (4.9) becomes 

^' + n^ + ^{n' -n^) = ^ip+p)a^(w-^y (4.24) 
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Eq. (4.10) becomes 



,A / . 9 , , (p 



_lxV + ix/-ia^(^-/^)=0. 



Eq. (4.11) (using the identity (3.62)) gives 



H' -H^ + -{p + p)a'^ 



+ 



X 



(4.25) 



(4.26) 



Finally, Eq. (4.16), using the background conservation equation (3.61) (with a = 0), can 
be rewritten as 



1 



1 



S' -{p + p)V^U - 3(p + :p) V + 3n{£ + P) - -^'{S - 3P) - -X'{p - 3p) = 0. (4.27) 



2'^ ' '2 

Eq. (4.17), using again the conservation equation, can be rewritten as 

1 



W' + HW - * + 



P + P 



p'W-P + l^'{p-3p){W-^) 

Z if 



(4.28) 



4.3 Relation to String-frame scalar perturbations 

In order to help the reader to make contact with previous results [11], it is useful to provide 
explicit relations between the E-frame and the S-framc perturbation variables. Unlike tensor 
fluctuations, scalar fluctuations are indeed neither gauge- nor frame-independent. Of course, 
since infinitesimal coordinate transformations have the same form in both frames, gauge- 
invariant perturbations can be defined in the same way in each frame, so that Eqs. (4.20), 
(4.21) can also be interpreted as the definition of gauge-invariant variables in the S-frame, 
provided the background and perturbations variables are referred to the S-frame metric, 
dilaton and fluid sources. 

The relation between a generic set of scalar perturbations in the E-frame (see for instance 
Eq. (4.19), completed by the fluid perturbations of Eq. (4.5)), and the corresponding set 
in the S-frame, can be easily obtained by recalling the transformation (3.53): 



QE = 

(Te) 



o,,e 



-<p/2 



HE = ns- ^'12. 



V 2ip 



(4.29) 
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We are considering the case c? = 3, and we have appended to our variables the labels E and 
s for an easier distinction of the two frames. It follows that 

Sg^. = Sg^.e-'^ - X Sgj^u, {5Te% = e'^{6Ts)l + 2x{5TEr^, (4.30) 
from which, using the definitions (4.19), (4.5) (valid in both frames), 

(t>E = (l>s - 2 ' ipE = tps+ 2 ' = Eg, Be = Bs, xe = Xs, 

6pE = Spse^"^ + 2xPE, 6pE = SpsS^"^ + 2xPE, ue = Ug. (4.31) 

For our uniform-dilaton gauge, in particular, xe = Xs = 0) so that all the metric perturba- 
tion variables coincide in the two frames, while the fluid perturbation variables only diflFer 
by an overall rescaling factor. 

For what concerns the gauge- invariant variables, using Eqs. (4.29), (4.31) we obtain, 
from the definitions (4.20) and (4.21): 

^E = ^s-\-^{B-E') = ^s-\Xs, Xe = Xs, 

^e = ^s + ^ + y{b-e') = ^s + \xs, nE = ns. (4.32) 

The (gauge-invariant) scalar field and curvature perturbations X and TZ are thus frame- 
independent, together with the (gauge-invariant) potential W, and the combination ^'-|-$. 
The two Bardeen potentials, on the contrary, are not separately frame-independent. 



5 Solutions for scalar perturbations 
5.1 General analytic considerations 

In this section we will compute the spectrum of scalar perturbations associated with the 
regular gravi-dilaton background (3.50). We will apply the equations obtained in the pre- 
vious section by setting everywhere p = p = 0, and = 5p = n = 0, as we shall consider a 
pure gravi-dilaton system with no fluid sources. 

We are interested, in particular, in studying the evolution and the final spectrum of the 
Bardeen potential * and of the curvature perturbation TZ. We thus need the decoupled 
equations describing the evolution of such variables. Let us first note, to this purpose, that 
a decoupled equation for ^ can be obtained by expressing X and its first derivative from 
Eq. (4.24), and using such result to eliminate X and X' from Eq. (4.25). The resulting 
equation can be written in the form 

^" + 2^ V + 2(rL' rL%\ - 4 ^^ ^ '^^^ = 0' (^-1) 
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where 



(5.2) 



It can be easily checked that, for a local potential, one has 4(7^^ — Ti.') = (p'"^ (see Eq. 
(3.57)), and one recovers from (5.1) the standard evolution of the Bardeen potential for a 
minimally coupled scalar field (see for instance [15]). In our case, the first derivative term 
in eq. (5.1) can be eliminated by introducing the variable U = and the equation can 
be rewritten in "pseudo-canonical" form as 

U" - ( clv' + 5 V = 0> (5-3) 



Z 

where, using the background identity (3.62), 



4 . . a2 dV 



cl = ^(W^ -n') = l + ^^e", 

l~i 'hi 
2a(W2 _ 7^/)i/2 = (^•'^) 

For what concerns the equation for curvature perturbations we may note that, in our 
gauge, 5(p = 0, so that TZ = —ip (from the definition (4.21)). The decoupled equation 
for TZ can thus be directly obtained from the decoupled equation for ip by combining Eqs. 
(4.12)-(4.15), and by using the background equations (3.58)-(3.60). The final result is 

z' 

n" + 2-7^' - c^v'^n = o, (5.5) 

where 

z = aif'/n, (5.6) 

and the first derivative term can be eliminated by introducing the variable v = zTZ, which 
satisfies the standard equation 

v" - (^clV^ + ^^v = 0, v = zTZ = -a(^X + . (5.7) 

It is important to stress that v, in the case of a local potential, exactly represents 
(modulo a sign, and in units Mp = 2), the gauge- invariant canonical variable diagonalizing 
the action for scalar perturbations [15]. In our case Eq. (5.7) is different from the canonical 
one because of the effective "speed coefficient" ^ 1. However, in our background (3.50), 

— 1 goes rapidly to zero at large negative times so that, even in our case, v is the 
appropriate variable for the asymptotic normalization of scalar perturbations to a vacuum 
fluctuation spectrum, namely 

\vk\ = (2A;)-V2 (5.8) 
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for 77 — > — oo. The normalization of Vk then fixes the asymptotic normalization of the TZj- 
modes, through the definition TZ = v/z. 

The normalization of TZ, in its turn, controls the quantum normalization of the Bardccn 
potential, once a relation between TZ^ and '^^ is given. Such a relation can be easily 
obtained again working in the uniform dilaton gauge, where TZ = —ip and = tp + TiE' . 
By eliminating ^ in Eq. (4.14) through Eq. (4.13), and Va^e'^ through the background 
equation (3.58), one is lead to the relation 

TZ' = = -^V'^r, (5.9) 

which holds exactly for all modes and at all times. For the asymptotic normalization at 
ri — —00, where all modes are oscillating outside the horizon {TZ^ ^ e"^''^), we then 

find, in particular. 

It follows that, for what concerns the /c-dcpcndcncc of the initial fluctuations, normalized 
to a vacuum fluctuation spectrum, one has TZk ~ f k^^^"^, while ^k ^ Uk ^ k^^^'^. 

It should be noted, finally, that in the case of a local potential = 1, so that Z = 
(see Eqs. (5.4), (5.6)), and the variables U and v satisfy the same canonical equation with 
"duality-related" pump fileds. Such a duality relation is broken, in general, for a non-local 
potential, but it becomes formally restored asymptotically, at large enough times. 



5.2 Analytic estimates using different matching prescriptions 

We shall now apply the above equations to the computation of the scalar perturbation 
spectrum associated with the background (3.50). 

We start noticing that for such background, sufficiently far from the bounce, the dilaton 
potential becomes negligible, and one recovers the usual vacuum solution of minimal pre-big 
bang models [1]: 

a ~ (— 7?)^/^, (/? ~ — A/31n(— 77), ?7 — s- — 00, (5-11) 
a-(r/)^/^ (p Vsin{r)), r) ^ +00, (5.12) 

It follows that, asymptotically, = 1, z"/z = a"/a, and Eq. (5.7) for v (or for TZ, 
Eq. (5.5)), exactly reduces to Eq. (4.4) for fi (or for h, Eq. (4.3)). We thus expect, at 
least for the low-energy modes exiting (and then re-entering) the horizon far enough from 
the bouncing region, that the curvature fluctuations TZk should be amplified with a scalar 
spectrum identical to that of tensor perturbations. 
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For such low-energy modes we know that the initial solutions, normalized to the quantum 
fluctuations of the vacuum, are given in terms of the Hankel functions of second kind and 
index zero, % ~ 7?.^ ~ (|r/^/^|/a)i7o d^'?!), which describe a logarithmic growth of the 
fluctuations outside the horizon. In the case of tensor modes, on the other hand, there 
seems to exist an almost unanimous consensus on the fact that the final spectrum should 
be computed by assTiming a smooth evolution of the metric fluctuation h across the bounce. 
In the background (5.11), (5.12) this leads to a low-energy spectrum which, after all modes 
re-enter the horizon (rj — > -|-oo), can be parametrized as 



k_ 

h 



Ai + A2 In 



(5.13) 



where ki is the high-energy cut-off scale, and ^1,^2 are constant (and /c-independent) 
dimensionless coefficients depending on the initial normalization. After the bounce, but 
still outside the horizon [kr] <C 1), the spectrum takes the form 



3/2 



[Bi + B2ln{kir])]. 



(5.14) 



For the scalar modes TZ^, of energy low enough to be unaffected by the specific kinematics 
of the bouncing region, we would thus expect a similar spectral behaviour, namely 



ki 



1 , k 
ai -I- a2 m — 

ki 



inside the horizon, and 



\5n\ = k'/^\nk\ 



ki 



\/3i+P2\n{kiri)] 



(5.15) 



(5.16) 



outside the horizon. Given S-ji, the corresponding Bardeen spectrum after the bounce is 
then determined by Eqs. (5.9), (5.10) outside and inside the horizon, respectively, as 



kri < 1, 



kri » 1. 



(5.17) 



The above expression for S-j^ and (5^ will be fully confirmed by a detailed numerical 
computation, as we shall see in the next subsection. Before presenting numerical results, 
however, we want to stress with an explicit example how the above spectrum is crucially 
dependent on the assumption that the variable TZ goes smoothly through the bounce, and 
how the alternative assumption of a smooth behaviour of ^' would lead instead to different 
spectra, in disagreement with subsequent numerical computations. 

Let us suppose, for our illustrative purpose, that in the absence of a full and explicit 
knowledge of the background dynamics near the bouncing transition we want to compute 
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the scalar spectrum by performing a matching of the perturbation variables, from rj = —rji 
to rj = +r]i, where rji is typically the transition time scale. Let us also suppose that, 
outside the transition regime (i.e. for rj < —rji and rj > rji), the contribution of the dilaton 
potential becomes negligible, and the background evolution can be correctly approximated 
by the aymptotic solutions (5.11), (5.12). We shall put, in particular. 



a = 



1 

ViJ 



\ 1/2 



H 

H 



1 

1 

2^' 



V 

V 



■■ -2V3a, 



v<-m, (5.18) 
v>+m- (5-19) 



We will concentrate our explicit example on the computation of the scalar spectrum out- 
side the horizon, and we will start our discussion by assuming a smooth behaviour of the 
curvature perturbation TZ across the bounce. In other words, the matching from —rji to r/i 
will be first performed by imposing the continuity of TZ and TZ'. 

We thus consider curvature perturbations of low enough momentum, outside the horizon, 
before the bouncing transition, and we solve (the Fourier transformed version of ) Eq. (5.7) 

for rjcx < rj < —rji, where {krjexl = 1- By using for the asymptotic expansion in the 
regime \krj\ <C 1 [18] (see also [1]) we obtain 



TZk{rj)=Ak 
z^TZ'kirj) = kBk 



l-k 



Vex 
2 



drj / z drj" 



Vex 



+ kBk 



Mjj' + • 



Vex 



r I" 

*' Vex Vez 







r' z-^drj" 




'Vex 





(5.20) 
(5.21) 



where A^, are integration constants to be determined by the initial conditions at ?7 = rjex- 
Such initial conditions, on the other hand, are prescribed by the canonical normalization 
(5.8) of the variable vj-, and the result is 



A. 



2k 



Bk 



2k 



(5.22) 



where a and /3 are appropriate phases, originating from (possibly) random initial conditions, 
and where Zex = — 2-\/3/ v^fc?7i (using the background (5.18)). The integration of Eqs. (5.20), 
(5.21) gives then, to lowest order. 



T^kiv) = Vrj^ [ci + C2 ln(-A:r?)] , 7^fe(^?) = C3v^i(-7?) S 



(5.23) 



where ci,C2,C3 are fc-independent complex numbers. It should be recalled that we are 
working in units M| = 2, and that the correct canonical dimensions of the Fourier modes, 
[7?.fe] = are easily restored when Eq. (5.23) is multiplied by \/2/Mp. 
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In order to obtain the spectrum of low-energy modes after the bounce, but before they 
re-enter inside the horizon, we wih now solve Eq. (5.7) for ryi < 77 < rjre, where krjre = 1, by 
using the values of TZ and TZ' computed dX rj = —rji, from eq. (5.23), as initial conditions 
at 77 = r/i. The asymptotic expansions (5.20), (5.21) are still valid, with the only difference 
that the lower limit of integration is now rji instead of rjex- The integration constants A^, 
are now determined by the conditions 

(5.24) 

(5.25) 



T^kim) =Ak= 7^jt(-77l) = [ci + C2 ln(A;77i)] , 

k 



T^'kim) = —Bk = T^'ki-Vi) = C3/V^ 
^1 

Using the background (5.19) we thus get, to lowest order, 
T^kiv) = VVi [ci + C2 ln{kr]i) + C3 ln(A;ir?)] , 



1 



12c3-^{krifln{km) 



Vl^V < lire, kfj < 1, (5.26) 



which provides the spectrum of curvature perturbations outside the horizon, after the 
bounce, in agreement with our expectations (5.16). The corresponding spectrum for the 
Bardeen potential is then obtained from Eq. (5.9), which leads to 



^kiv) 



{kriY 



2'' 



{krjf 
16 



ln{kr)i 



(5.27) 



The leading term, for krj <C 1, is in agreement with the expectation (5.17). 



The above results for the scalar perturbation spectrum, obtained through a matching 
procedure based on the continuity of TZ and TZ', will be confirmed by the numerical compu- 
tations of the next subsection. Let us now show that a matching based on the continuity 
of * and leads to different results for the scalar spectrum. 

Following the same procedure as before, we start considering the low-energy modes of 
the Bardeen potential outside the horizon, before the transition, and we solve Eq. (5.3) for 
Vex < ?y < ~Vi- By using the background (5.18) we are lead to the asymptotic expansion 



1-k^ Z-^dv' / Z^dr]" + 



df]' r z^df]"] + kBk r 

Vex I fjex 

k'^T] 



Z-^dr]' 



, (5.28) 



Vex 



Vex 



2^2 



Z^dv' 



Vex 



(5.29) 



where A}-,B}- are determined by the normalization of Uj^ at 77 
\iex^k{riex)\ = k-^l"^. The result is 



A. 



2^/3 
k'^VVi 



Bk 



kV3 V 



rjex, namely \Uk\ = 
(5.30) 
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We thus obtain, to lowest order, 



(5.31) 



{k7]y ' "^ '^ k^r]3 ' 

where 61,62 are fc-independent complex numbers. It may be noted that, before the bounce, 
the spectral distribution of is in agreement with that we would obtain from the spectrum 
of TZk, eq. (5.23), by using the general relation (5.17). 

Let us now compute 'ffc after the bounce, by solving eq. (5.3) for rji < i] < rjre before 
re-entering, and using ^'fc(— r?i) and ^^,(— r/i) as initial conditions at r] = rji. The asymptotic 
expansion after the bounce, using the background (5.19), and taking into account all the 
required flips of sign, takes the form: 



m 
4^2 



Ak(l-k^ r Z-^dri' r Z^dr]") + kBk P Z'^dr]' 

i-e r z-'dr,' r' z'dv" - r 



(5.32) 



Z'^dri 



m 

'2rf 



Bh 



rz-'dv' 

Jrii 



+ 



krj 
2Z2 



1-k' I Z'^dr]' I Z-'^dri' 



Vl 



-2 



(5.33) 



The integration constants are determined by the continuity of ^' and ^' as 



3/2 



2 ' 



1 3/2 



(5.34) 



The super-horizon spectrum of the Bardeen potential, after the bounce, to leading order, 
is then determined by 



*fc(r/) 



{k'qY 



hi 



bn 



{kmy 



r] > r]i, kr] <C 1. 



(5.35) 



This shows that the Bardeen modes, in case of a smooth crossing of the bouncing region, 
would tend to stay constant after the bounce (outside the horizon), instead of decreasing 
in time like rj^^, as predicted by Eq. (5.27). 

Let us now compare the above results with a numerical integration of the perturbation 
equations. 



5.3 Numerical analysis and comparison with analytic estimates 

For a numerical computation of the spectrum it is convenient to express the evolution of 
our perturbations TZ and ^ directly in the S-frame, where we know the explicit form of the 
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background solution, Eq. (3.50). In terms of S-frame variables, using the gauge Sip = 0, we 
have, in particular (see Subsection 5.2), 



TZe = -ipE = -tps, 

^E = ^E + UeE'e = ^s+ {Us - K. (5.36) 

^From the system of scalar perturbation equations (4.12)-(4.15) we can then obtain a system 
of first order, coupled equations for the variables ■i/'s and i?', whose numerical solutions will 
allow us to determine TIe and ^e everywhere. We will find that 7?., unlike goes smoothly 
through the bouncing transition, and that the numerical results provide a good fit of the 
scalar spectra analytically obtained using the continuity of Tl and Tl' . 

Dropping, for simplicity, the subscript referring to S-frame variables, and using the 
transformations (4.29), (4.31), the perturbation equations (4.8)-(4.10) can be written in 
the S-frame as 

E" + {2n - ip')E' + - </) = 0, (5.37) 
2V'' + (2H - (/?')</' = 0, (5.38) 
2VV + (2H - ^')V^E' - 3V''(2W - if') - {m^ - m(p' + cp'^) = 0. (5.39) 

By eliminating cf) through Eq. (5.38), and moving to the S-frame cosmic time coordinate 
t, the combination of the above equations leads to the first-order autonomous system of 
equations: 

i];k = Ak{t)i>k + Bk{t)rk, (5.40) 
tk = Ckit)iJk + Dkit)Tk, (5.41) 

for the variables V'fc = —T^k a^id = a^E^ (an ovcrdot denotes differentiation with respect 
to t). The coefficients appearing in Eqs. (5.40) and (5.41) are 

2(y-2g) 2 f ^-2H Y 2 

^ \ (p J 

Ck{t) = ^-h Dk{t) = ip-H-Ak{t). (5.42) 

where uj = k/a denotes the physical momentum. Using eq. (5.36) it is also clear that the 
E-frame Bardeen potential is given by 

*f = V'fe-^(</'-2iJ)Ffc. (5.43) 

In order to solve numerically Eqs. (5.40) and (5.41) in the case of the non-singular 
model introduced in the previous sections, quantum mechanical initial conditions are to be 
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imposed at early enough times, when the relevant mode is still inside the horizon. We define 



7^fe = Ffe + iQk, nk = Gk + iPk, (5.44) 

w{t) = VEz = (5.45) 

and we use the canonical normalization (5.8). In particular, we choose the initial phase at 
t = ti in such a way that the initial value of TZk is real, namely 

FkiU) = )=, QkiU) = 0, 

Gk{U) = (^-t-) -A^, mU) = -M^. (5.46) 



Using (5.40), (5.41) the initial conditions for the real and imaginary parts of are then 
fixed as follows: 



Mrfe(*i)} = - ( 



Gk - AkFk \ 



Im{r,(t,)} = - ('^M^l • (5-47) 



-A 

Bk /ti 

The spectrum of all the interesting quantities can now be obtained through a numer- 
ical integration of Eqs. (5.40), (5.41). Consider, for instance, the spectrum of curvature 
perturbations S-ji = k^/'^\Rk\- In Fig. 1 we have illustrated our numerical results for the 
super-horizon evolution of S-r, for different comoving momenta, normalized to a quantum 
spectrum of vacuum fluctuations. We have used the background solution (3.50) with to = l, 
and we have fixed the absolute amplitude of the normalized spectrum by choosing <^o in 
such a way that 

{^) — = 10-- (5.48) 

The mode n = kto ~ 1 starts its first oscillation already for t ~ 0, so that it basically exits 
the horizon and then immediately re-enters. Smaller wave-numbers stay longer outside 
the horizon. We have only considered low-fequency modes, k <^ 1, which are consistently 
amplified in the given model of bouncing background. Higher frequency modes may develop 
a quantum gravitational instability (near the origin t = 0, where < 0), similar to that of 
tensor modes in higher-derivative string cosmology backgrounds [19]. For such modes our 
(low-energy) model of bakground is clearly inadequate. 

^From Fig. 1 it is already apparent that the spectrum of TZ, outside the horizon, follows 
the predicted power-law distribution k^/^, in agreement with Eqs. (5.16) and (5.26). For 
a more accurate determination of the spectrum we will choose a fiducial time = 100, 
for which all modes with k < 10~^ are still outside the horizon. The spectrum |(57e(Ta;)| 
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Figure 1: Numerical results for the time-evolution of ^7^, and for different values of the 
comoving wave-number (we have set 1^ = 1, and used the normalization (5.48)). 

can then be computed as a function of k, and the obtained result can be fitted with the 
analj4;ical expression of STZ}- obtained with the matching procedure outlined in the previous 
subsection. 

For numerical reasons, it is convenient to work with the squared modulus of Eq. (5.26) 
(since the coefEcients appearing in that equation are complex). We thus set, at r = r^;, 

\5n(T:c)? = K^(ai + A luAt + 71 In^ k), (5.49) 

where ai, (5\ and 71 are real numbers. In Fig. 2 the numerical results for the values of 
are reported with the open diamonds, for a selection of modes 10~^ < < 10~^, and 
Tx = 100. Such values are well fitted by Eq. (5.49) with 

ai = 0.339 X lO-'', /3i = -0.345 x 10"^, 71 = 0.664 x 10"^ (5.50) 

The corresponding curve (plotted in Fig. 2 with a dashed line) nicely intepolates through 
the numerical points, confirming the previous results for the analytical estimate of the 
spectrum of 7?.^. Note that, when the spectral amplitude is expressed in terms of the 
correct dimensional factors, we obtain (in our model of matching) ai ~ /3i ~ (iJi/Mp)^, 
where Hi = H(r]i). Eq. (5.50) thus represents a fit of the matching scale Hi, in agreement 
with the assumed normalization (5.48). 

A similar numerical analysis can now be repeated for the spectrum of the E-frame 
Bardeen potential, from Eq. (5.43). The numerical solutions show, first of all, that 
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logK ^ 



Figure 2: Spectral distribution of \Stz\, at the fixed fiducial time = 100 when all the 
selection of modes is well outside the horizon. The dashed line is the fit given in Eqs. (5.49) 
and (5.50) (we are using decimal logarithms on both axes). 

the values of ^ and ^' (though continuous at t = 0) change drastically while crossing the 
bouncing region (unlike TZ), as illustrated in Fig. 3. This support the expectation that the 
correct spectral distribution shoTild be obtained, analytically, by matching the curvature 
perturbation variable, instead of the Bardccn potential. 

This expectation is indeed confirmed by Fig. 4, where we plot the Bardccn spectrum, 
as a function of time, for different comoving momenta which are all inside the horizon. 
The "post-bounce" evolution of the given super-horizon modes, instead of being time- 
independent as required by Eq. (5.35), decreases in time exactly as predicted by Eq. (5.27), 
namely '^kiv) ~ 'n~^- Sufficiently far from the bounce, where the S-frame cosmic time t is 
related to the conformal time ij by 

t = J asdr] = J aEe'^/'^dr] ~ r]^^+^y^ (5.51) 
(we have used eq. (5.12)), we thus obtain 

^,y-2 ^r^/(3+^/3)^ 4/(3 + x/3) ~ 0.845. (5.52) 

This general behaviour predicted by Eq. (5.27) is represented in Fig. 4 by the dashed 
curve, which nicely agrees with the time-evolution of the Bardeen spectrum for all plotted 
(super-horizon) modes. 



30 



Figure 3: Evolution of the (E-frame) Bardeen spectrum around the bounce. 




Figure 4: Evolution of the (E-frame) Bardeen spectrum for super-horizon modes, after 
the bounce. The dashed curve parametrizes (with an arbitrary normalization) the time 
behaviour rf'^ ~ ^-4/(3+73) ^ ^-0.845 predicted by Eq. (5.27). 
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Figure 5: Spectral distribution of |(5^| at the fixed fiducial time Tx = 100, when all the 
modes included in the plot are well outside the horizon. The dashed line is the fit reported 
inEq. (5.53). 

Also in agreement with the analytical prediction (5.27) is the spectral distribution of the 
super-horizon modes reported in Fig. 4 , which can be easily inferred from the (decimal) 
logarithmic scale of the figure as a more precise computation 

of the spectrum we may observe that the numerical values of \S^\, at a given fiducial time 
Tx = 100 such that all the considered modes are still outside the horizon, can be fitted by 

= -T72 , «2 = 0.347 X 10"^ (5.53) 

Such a numerical fit, illustrated in Fig. 5, precisely confirms the leading order distribution 
analytically predicted in Eq. (5.27). 



6 Conclusion 



If taken at face value the results presented in this paper indicate that the curvature per- 
turbation 7?., rather than the Bardeen potential behaves smoothly through the bounce. 
The same conclusion can be reached by looking at differential equations obeyed by ^ and 
7?., Eqs. (5.1) and (5.5), and at their formal solutions Eqs. (5.20) and (5.28), respectively. 
The solutions for H, and its canonically conjugate momentum I!-;^ = z^Ti! only contain the 
backgrounds (through z) under integration, thus smoothing out their rapid variation at the 
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bounce. By contrast, the solution for * and its time derivative contain, besides integrals of 
and additional rapidly varying pre-factors. The only way to make ^ smooth across 
the bound appears to be through the insertion of (5-function like contributions to Z so that 
the variations of the prefactors are carefully cancelled by the integrals. Such (5-function 
contributions will instead make 7?. behave wildly through the bounce. 

In particular, the fact that I!-/^ and * cannot be simultaneously smooth across the 
bounce stems from the simple relation between them 

k^^ - Ha-^Un , (6.1) 

(see Eq. (5.9)), and the fact that the quantity Ha"'^ clearly changes sign across the bounce 
(recall that we are using Einstein-frame quantities). If, on the basis of our example and 
of the above discussion, one concludes that TZ is continuous across the bounce, the spectra 
of adiabatic scalar perturbations in bouncing cosmologies is expected to be far from scale- 
invariant. In this case the only known way to rescue the phenomenology of PBB cosmologies 
is through use of the so-called "curvaton" mechanism [20] (see for instance [1, 21]). 

On the other hand, one has to concede the fact that our mechanism for producing the 
bounce through a non-local potential tends to tie its occurrence to a critical value of the 
scalar field, and thus to privilege continuity of perturbations on comoving (constant (p) 
hypersurfaces. The proponents of a scale-invariant spectrum [4] claim that the situation 
is totally different in the ckpyr otic/ cyclic scenario, provided this is looked at from a five- 
dimensional perspective. Possibly, such an approach will fix the matching conditions in a 
completely different way: what appears to be fine-tuning in four-dimensions could appear 
to be natural from a higher dimensional point of view. As a result, even the spectrum of 
tensor perturbations claimed in [13] differs from the one proposed earlier in [3, 22] by the 
lack of a logarithmic behaviour of the Bogolubov coefficient. 

In conclusion, we believe it would be very interesting to complement the work of [5] by 
linking in a precise way different matching prescriptions to different dynamical mechanisms 
for inducing the bounce and avoid the singularity. This will eventually connect physical 
observables, such as the tilt in the power spectra, to the very non-trivial, and yet poorly 
understood, dynamics that supposedly allows one to evade the Hawking-Penrose singularity 
theorems and thus to conceive the possibility of a pre-big bang Universe. 
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